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Abstract 

By using the solution describing a black hole embedded in the 
FLRW universe, we obtain the evolving equation of the black hole 
mass expressed in terms of the cosmological parameters. The evolving 
equation indicates that in the phantom dark energy universe the black 
hole mass becomes zero before the Big Rip is reached. 



1 Introduction 

In recent years, strong observational evidence shows that the current ex- 
pansion of our universe is accelerating [I]. In the framework of the general 
relativity, this positive acceleration means that, at present, our universe is 
dominated by a mysterious component, the dark energy [2]. Usually, we may 
model the dark energy by a perfect fluid with a negative pressure p = wp < 0. 
The dark energy may be in the form of a vacuum energy with w = 1, or a 
quintessence model with w > —1. The observation gives the constraint [3]: 
w = — 1.07lQ;o9(stat l(T)±0.13(sys). Then there exists the possibility of the 
phantom models with w < —1. In this case, the universe undergoes the 
"supper accelerated" expansion. Finally, in the finite time, the Big Rip sin- 
gularity is reached, and all bounded objects in the universe are torn apart 
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What would happen to black holes in an expanding universe filled with 
the phantom dark energy? In [S], E], using the Schwarzschild metric, the 
authors have obtained the result 

M ^ ^ = At^G^AM^Ip^ + Poo], (1) 

where poo and poo are respectively the pressure and energy density of the 
universe at the asymptotic limit r +00. And A is a dimensionless constant. 
Then the author concluded that the masses of all black holes will tend to 
vanish as the universe approached the Big Rip. In [7], the authors obtained 
the same result ([T]), by using the non-static Schwarzschild metric 

_ (1 _ 2G^^j ^ _ 2GK^y ^ ^ 

(2) 

However, the metric is asymptotically fiat and does not describe exactly the 
space-time of a black hole embedded in the FLRW universe. Actually, Eq. ([T]) 
is obtained by ignoring the backreaction of the phantom matter on the black 
hole. When the matter density of the background is low, the ignorance 
is reasonable. However, as the Big Rip is approached, the matter density 
becomes very large. The backreaction should not be neglected, and the 
metric ^ becomes unsuitable for describing the spacetime. Then, we should 
adopt some exact metric, in order to take into account of the backreaction. 

Recently, new exact solutions describing black holes embedded in the 
FLRW universe have been suggested in [HI [9] . Then it is interesting for us 
to study the problem of the dark energy accretion on black holes by using 
the exact solutions. In [S], the accretion rate (Hawking-Hayward quasi-local 
mass added per unit time) is obtained. But the result is not suitable for the 
application in the cosmology. In order to study the evolution of the black hole 
mass as the universe expands, we should express the accretion rate in terms 
of the parameters at the asymptotic limit r — -|-cx), namely the cosmological 
parameters, Poo, Poo, et al.. 

In this paper we will try to obtain such an expression by using the solution 
with imperfect fiuid and radial mass fiow in [8J. Our result shows that the 
evolution of the black hole mass is governed by a second-order differential 
equation, which is different from Eq.([T]). 

The paper is organized as follows. In Section [21 we recall the results in [S]. 
In Section [31 we show the evolving equation of the black hole mass expressed 
in terms of the cosmological parameters. Section [H contains a discussion and 
conclusions. 
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2 the Solution with Imperfect Fluid and Ra- 
dial Mass Flow 



In [8], the line element is written in the form 



(3) 



1 — ^^^7^- Here and after, we take 



A^{t,r) 

with A{t,r) = 1 + ^ and B{t 
c = h = 1. Then, at the asymptotic limit r +oo, the metric becomes 
the FLRW metric. In this space-time, the physically relevant mass is the 
Hawking- Hay ward quasi-local mass [HI [IDl flT] 



'i^uif) = a(t)m(t). 



(4) 



In this paper, we study the solution with radial mass flow, u{r,t). The 
cosmological matter is assumed to be described by the imperfect fluid energy- 
momentum tensor 



T"'' = {p + p)u''u' + p^"" + g'u" + q% 



ab I „a„,b 



J>r. A 



with 



= (^^l + aM%^ M, 0, 0^ , g'^ = (0, g, 0, 0). 



(5) 



(6) 



Here, the purely spatial vector field q"- describes a radial energy flow. The 
condition 

lim u{t,r) = lim q{t,r) = 0, 

r— >+oo r— >+oo 

should be imposed, in order for the momentum-energy tensor to become 
the perfect fluid momentum-energy tensor in cosmology at the asymptotic 
infinity r ^ +oo, 

T£ = (Poo + Poo)w^M^ + p^g^^. 
The Einstein equations give four independent equations [S], 



rriH 
SttGp 

8nGp 
Q 



--aB^AvT+oMHi^ip + p)u, 



(7) 



3C^ + 2{C + 



GmC. 



+ {C + 



rAB ' _ 
GfnC^ la^A^u^ 



rAB ' 1 + aM%2 



-b + p) 



u 



(9) 
(10) 



with A = j J d6dipy/g^ = Anr'^a'^A'^ representing the area of a spherical 

surface with isotropic radius r, and C = f + Here and after, the dot 
denotes the derivative with respect to the time. 

This is the solution describing a black hole embedded in the FLRW uni- 
verse [8]. As r oo, Eq.Q becomes the Friedmann equation 



.a' 



a2 



Then Eq.(l9]) may be consider as a generalized Friedmann equation. Eq.Q 
and ([9]) indicate a generalized conservation law 

p(l + a^AV) + ^(a^AV) + 3(- + 24) (p + p) (1 + a'^V) 

a A 

+ (j9 + p)2a2AV(- + 24 + -) = 0, 
a A u 

which, as r — i> +00, becomes the conservation law in the FLRW universe 

p + 3-(p + p) = 0. 
a 

Additionally, the radial component of the energy- momentum conservation 
law, Ti|^ = 0, gives another conservation law 

= J 5(^ + 2^){l + aW)(p + p)u + ip + p)a^AV(^ + 4 + ' 

2B\J\ + a2y4%2 a A a A u 

+ (1 + a^A\^) {{j> + p)u + (p + p)u\} + p' + - + + a'^V) 

(12) 

Here, we note that the time component of the energy-momentum conserva- 
tion law, Tof^ = 0, does not give an independent equation, which can be 
derived by using the equations ([7]) and (fTTj) . 



3 the Evolving Equation in terms of Cosmo- 
logical Parameters 

Eq.© gives the accretion rate. Since -u < 0, the mass mii increases if 
p -|- p > 0, stays constant in a de-Sitter background, and decreases for the 
phantom dark energy p < 0. However, this equation is not convenient for 
us to study the evolution of the black hole mass as the universe expands. In 
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cosmology, what we need is an evolving equation of mn expressed in terms 
of the parameters at the asymptotic limit r — > +00, namely the cosmological 
parameters. 

To obtain such an evolving equation, we may take the limit of Eq.(I71) as 
r ^ +00. Then we have 

rriH = -^TTa^lPoo + Poo) hm (wr^), (13) 

where Poo(^) and poo(^) are respectively the limit of p{t,r) and p{t,r) as 
r ^ +00. This result indicates that, as r — > +00, 

«(t,r) ~ ^^^ + 0(r-3). (14) 



r 



Similarly, as r — » 00, we may expect 



p(t,r)^p„o(t) + ^ + 0(r-^), (15) 
r 

p(t,r)^Poo(t) + ^ + 0(r-^). (16) 
r 

Then, up to the order r~^, Eq. ffT2|) gives 

^[a^(Poo + poo)uj\ = 2a^[pi - Gm{poo + Poo)] (17) 
At the same time, up to the order r~^, Eq.l fTTl) indicates two equations 

Poo + 3-(poo + Poo) = 0, (18) 

ct 

pi + 3-{pi + Pi) + 3Gm{p^ + p^) = 0. (19) 
a 

Here, naturally, we can impose the equation of state 

Pooit) = wpooit), (20) 

where w is a dimensionless constant. Further, we assume the same equation 
of state for pi{t) and pi{t) 

pi{t)=wpi{t). (21) 
Together the two equations and Eq. (|TSl) . we can solve Eq. (|T^ . and obtain 

pi(t) = 3^(7710 - m)(poo + Poo), (22) 



where mo is an integral constant with the dimension of mass. Now, Eq. ffT7|) 
can be rewritten as 

^[a^(Poo + Poo)moo] = 2Ga^[3wmo - (3w + l)m]{j>^ + p^). (23) 

Using this equation and Eq. (fT3|) . finally we can obtain the evolving equation 
of ruH 

ffiH + '^~^H - 4:7tG[{3w + l)mH - 3wmoa]{poo + Poo) = 0, (24) 

where rriH = ^ • This is a second-order differential equation as we expect 
from the second-order differential equations ([7])-([ni)- 



4 Discussion and Conclusion 

In the last section, we obtain the evolving equation of the mass of a black hole 
embedded in the FLRW universe. Both Eq.([l]) and our result, Eq. (1241) . imply 
that the black hole mass decrease in the phantom dark energy universe, but 
the qualitatively evolving behaviors of the black hole mass indicated by the 
two equation are different. Due to Eq.([T]), it has been shown in [31 [12] that 
the mass of a black hole tends to zero in the phantom dark energy universe 
approaching the Big Rip. However, by using our result Eq. fl^ . we find that 
the black hole mass in the dark energy universe has been zero before the Big 
Rip singularity is reached. 

In order to show this, firstly let us rewrite Eq. in terms of the deriva- 
tive with respect to the scale factor a(t), 

d^ruH , 3(1 -w) _^^dmH 3 , in/i , ^ -2 ^ , - _i 

, ^ H a — — = -[3w + + w)a mn - -w[l + wjuiQa . 

da'' 2 da 2 2 



(25) 



Here, we have used the equations 



'^^^ —Poo, Poo oc a-^(^+-). 



The general solution of Eq. ( l25i) is 

mnit) = CMt)]''-''" - C2[a{t)]-''-'+^^/' + ^l±^moa(t), (26) 

where Ci and C2 are the integral constants. Use to to denote the initial 
moment. The two constants are determined by the initial mass, mH{to), and 
the initial accretion rate mH{to)- 
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In a phantom dark energy scenario, w is close to —1 and w < —1. Then 
we may have 

3(1 + w) 

l + 3w<0, l>-31 + w/2>0, ^ ^<0. 

3w + 5 

Using muito) > and muitQ) < 0, generally we may have 

Ci > 0, C2 > 0. 

Then Eg. (1261) indicates that, before the Big Rip singularity is reached, the 
black mass has been zero at the moment te which is determined approxi- 
mately by the equation 

Physically, this result seems to be better, since the Big Rip singularity be- 
comes unnecessary for the disappearance of the black holes. Of course, the 
analysis above may be not rigid when the black hole mass is at the order of 
the Planck mass, since the effect of the Hawking radiation might dominate 
the evolution of the black hole with the Planck mass. 

In the author have shown that Eq.([l]) can be generalized to the non- 
static metric ([2]). Then we might expect that Eq.([T]) should be reduced by 
taking a(t) to be constant, d = 0, in Eq. ip^ . However, it is not the case. 
The reason is that Eq.([l]) in [7] is obtained under the assumption 

lim (ur^) = AG^M^. 

r— >oo 

From Eq. (JT711 . we know the assumption for the solution used in this paper is 
unreasonable. Then it becomes natural that Eq. ([T]) can not be reduced from 
Eq.d^ij) with d = 0. 

Even, we should note that the previously known solution for the perfect 
fluid, Eq.([T]), cannot be reduced from our result Eq.([21]) by taking q'^ = 
in Eq.(l5]). In fact, g'* 7^ is necessary. It has been shown in [8] that, for 
the perfect fluid described by taking = in Eq.(l5]), the Einstein equations 
indicate 

p + p = 0. 

This can be also learned from Eq. ffTUj) . This implies that only the de-Sitter 
equation of state p = —p is allowed. Eq.(I7]) accordingly yields rfiH = 0. So 
no accretion happens for a single perfect fluid. Then the imperfect fluid ([5]) 
with g'^ 7^ is necessary for the accretion to happen when the metric ([3]) is 
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used. This indicates that the previously known accretion rate ([T]) is not the 
hmit of our result ([211) as — > 0. 

However, We do not think our result (l24l) is better than the result ([T]). The 
reason is that the imperfect fluid solution which our result ( 12^ is based on 
has the unphysical behavior [8] — the superluminal motion of the fluid as r — 
m/2. For this reasonour result cannot be used as a "better approximation" 
as compared with the result ([T]). We may take Eq. (!^ as an interesting try. 

Summarily, in the paper, by using the exact solution describing a black 
hole embedded in the FLRW universe [8] and the equation of state ( |20l) . 
we obtain the evolving equation of the black hole mass expressed in terms 
of the cosmological parameters, which is convenient for the application in 
cosmology. Our result, a second-order differential equation ([211), is different 
from the result in [Kj. Our result indicates that the black mass in the 
phantom dark energy universe has been zero at the moment te before the 
Big Rip moment ts, while, due to Eq.([l]), the mass of a black hole tends to 
zero in the phantom dark energy universe approaching the Big Rip. Here, we 
emphasize again that our result is not better than the result ([11) because 
of the unphysical behavior of the imperfect fluid solution which our result 
(l24l) is based on. 
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